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ABSTRACT 


An empirical study of pharmacy item demand from the 
Supply Department in a Naval Hospital was conducted to test 
three stochastic inventory models to see if investment in 
inventory could be reduced. The three models, A3RSR Method 
(Tukey) , Exponential Smoothing Method, and Regression Fore- 
casting, were run with eighteen months of demand figures to 
develop the forecasts and compare the models. The results 
were evaluated on the basis of numerical differences and 
dollar values of the forecasts vs. the actual demands. 

Exponential Smoothing proved to be slightly better, but 
no great differences were noted between the performances of 
the three models. One possible explanation was found to be 
excessive variance in the demand generated by the Pharmacy, 
due to inputs other than actual customer demands.  Recom- 
mendations were made for the institution of a basic system 


to aid the Pharmacy in determining their actual demand. 
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ROD C LUN 


The ability to forecast demand is a vital part of any 
inventory management system. Generally such systems answer 
two questions about the items they manage; 1) How much to 
order, and 2) When to Order. {1 Forecasting assists in answer- 
ing the first question (How much to order), therefore the 
ability to forecast accurately has a direct bearing on 
the performance of the total system. 

Forecasting, as defined by Brown, involves the "projection 
of the past into the future", while prediction relates to 
iu cCiDation of changes and new factors affecting demand". 
This separation of forecasting and prediction is important 
in that it binds forecasting to patterns of demand that have 
occured in the past. Thus, in Brown's view, forecasting has 
a historical background, whereas prediction has no real 
basis other than judgement and opinion. Among other experts 
in the field of forecasting, there is no absolute agreement 
AREAS point, and in fact the dividing line is at best 
obscure. Forecasting and prediction may be likened to 
objective and subjective methods of future demand determina- 
tion; the subjective methods employ processes which are not 
well specified and not easily replicable, while objective 
methods use well specified processes, can be computerized. and 
can be easily redone by others if necessary. ? For purposes 
of inventory management in a medical enviroment, the historical- 


based forecasting methods are generally preferable to 





subjective predictions. 

Within medical facilities, proper operation of the supply 
function is a necessity; stockouts could affect the mission 
Of the activity and impair medical treatments. To guard 
against stockouts, the manager may elect to hold a certain 
level of items, commonly called a safety level, in reserve 
against the possibility of a stockout. The safety level is 
intended to provide a supplemental supply of items in excess 
of the forecast if the demand should exceed that forecast. 
The safety level increases the costs associated with holding 
that inventory, but as long as this cost is less than the cost 
Of a stockout, the manager will elect to incur it. In this 
context, forecasting becomes doubly important. The medical 
inventory manager must be especially diligent about how he 
determines the amount of material he buys and stocks. The 
penalty for a stockout may be extreme, yet the manager 1s 
constrained from stocking material at a prohibitively high 
safety level by monetary, space, and manpower limitations. 

He must find a forecasting system that allows him to stock 
material at an acceptable safety level, yet minimize the costs 
associated with holding the material. 

The rational supply manager will recognize the possibil- 
ities that a forecast may be either above the real demand 
observed, exactly equal to the demand, or fall short of the 

ctual demand. Within stock fund limitations imposed by the 


Fleet Material Support Office, if he is above or exactly equal 


to the demand, he need not fear a stockout. (Navy medicai 


supply managers finance their inventories through a working 








capital fund called the Navy Stock Fund, which is administered 
and authorized on an activity-by-activity basis by the Fleet 
Material Support Office. This Office analyzes activity stock 
fund requirements quarterly, and grants Stock Fund monies 
based on these analyses. Each supply manager receives this 
analysis, the Regular/Supplemental Grant Worksheet, with the 
quarterly Stock Fund Allotment/Suballotment Authorization, 


and can ascertain the reasoning behind the stock fund grant) ^ 


This paper will explore various forecasting models and 
their performance in stochastic demand situations. The demand 
data were gathered from the Naval Regional Medical Center, 
Oakland, California, and represents demand over time for items 
used by the Pharmacy in Federal Stock Class (FSC) 6505. These 
items were selected because they are among the most heavily 
used items in a medical facility, and are of paramount impor- 
tance in patient treatment. As stated earlier, this paper's 
scope is limited to forecasting models, and will not address 
other inventory related topics such as reorder point deter- 
mination, lead times, inventory cost determinations, and 
Economic Order Quantities. The reason for this limitation 
RE Cope is that forecasts are vital requirements of all 
inventory management systems; both civilian and military 
applications need them to carry out their functions. Also, 
forecasts can be tested empirically, with the output being 
applicable to that particular system. 

In the military medical supply environment, the systems 
in use are generally simpler than other systems in use through- 
out the rest of the Supply System. The problems faced by the 
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Supply System in general are more complex (repairable items, 
sizeable numbers of reservation items, etc.) than the medical 
supply problems, and therefore the models at higher System 
levels are more sophisticated. The models in use at the 
medical activity level tend to be simpler, general models 

Of the type used in managing smaller inventories. For instance, 
the model in use in medical activities without sophisticated 
data processing capabilities (manual labor intensive systems) 


is given by: 


l 


EEE 2 = Reorder Point 


Annual Demand Quantity x 5 = Requisitioning Objectives? 


This model is representative of the type used in a wide variety 
of civilian and military inventory systems for manual and simple 
E.A.M. applications. The first part of each equation, Annual 
Demand Quantity/12, is a moving average that gives a forecast 
of demand for the next month based on the monthly average for 
the past twelve months. The second part of each of the two 
equations, multiplications by two and five respectively, give 
inventory maintenance levels based on that forecast. 

It can be seen that the forecast is an element of inventory 
management that is essential to all systems; analysis of an 
inventory topic not common to all systems would have limited 
the study's usefulness. It was for this reason that the 


decision was made to limit the study to forecasting. 





PROBE DEE EN ET LON 


At the beginning of the study (and currently), the fore- 
casting model in use at the Naval Regional Medical Center, 
Oakland, California was the following; 


CD + 2(OP,) 


= OP2 
3 


CD = Current Demand 
OPı= Order Point Last Month 
6B>= Order Point Current Month 
Many questions can be asked about this model; 
1) Who had originally formulated and implemented the 
model? 
2) Had any studies been done prior to implementation? 
3) Would another stochastic model do better? 
Interviews with the Supply Officer of the Regional Medical 
Center revealed that the origin of the model was completely 
unknown. No one in the Supply Department or the Data Processing 
Division knew who had originally programmed the model. The 
model was a part oí the reorder program run for the Supply 
Department by the Data Processing Division, and had been in 
use for a long period of time, with even its original imple- 
mentation date unknown. The Supply Officer, in an effort to 
learn more about the model and its effects on costs and stock 
levels, had to request a print out of the reorder program, 
then analyze it in depth to understand the basic format of the 
model set forth above. It was quite evident that the fore- 


casting model in use was poorly understood in terms of its 


10 





origin and effects on supply operations. 

The effects on costs and stock levels of the model forecasts 
were also unclear. With no other data to evaluate then the 
actual performance of the model, no comparisons could be made 
and no judgements arrived at. The Supply Officer, in an effort 
to lower costs and improve performance relative to supply 
levels within his department, sought more information concern- 
ing the model he was using. A study of the model in use, along 
with the modification to it and other stochastic models that 
could be used, was decided upon as the best vehicle to see if 


improvements were possible. 


LL ANALYTICAL APPROACH 


To address the problems of how well the system performs 
mn relation to other stochastic models, what the effect of 
other changes to the present model would be, and how well the 
present and other models do in meeting the demand, an empir- 
ical study of the present system and alternate forecasting 
models was undertaken. An empirical study was decided upon 
primarily because it would use the inputs (demand, unit prices) 
of the actual system to be evaluated and reflect the perform- 
ance of forecasting models relative to the actual system. 

The framework of the study included three steps; 

1) The selection of a sample of inventory items to use 

in conjunction with forecasting models under study 

2) An analysis of forecasting models in actual operation 


with the sample 





3) An evaluation of the forecasting model performance 

in a managerial context 

The items carried by the Pharmacy at the Naval Regional 
Medical Center were selected as the source of the sample, as 
they represent some of the most important and heavily used 
items in the hospital setting. FSC 6505 items (drugs) are 
held in stock by the Supply Officer, and are drawn by the 
Pharmacy for the filling of prescriptions. Within the 
hospital setting, the Supply Officer represents the system 
supply of items available upon demand by the Pharmacy. The 
Pharmacy, as the end-user of the items, creates the demand for 
the FSC 6505 items when it requisitions them from the stock 
the Supply Officer holds. 

The system to be empirically tested sets the reorder 
amounts in the Supply Department to insure material avail- 
ability for the Pharmacy and other supply users. The Pharmacy 
generates demand through the requisitioning process, anda 
history of these demand figures are retained and used by the 
Supply Officer. The forecasting model uses these demand 
figures to forecast demand in the upcoming period (one month), 
and procurements of material for the Supply Department are 
made based on the forecasts. The value of the forecast is 
measured by how well the forecast foretells the demand and 
prompts procurement action to meet it. 

A demand history of twelve months accumulated demand recorded 
in the Retail Asset Stock Card (RASC) was used as the basis 
for sample selection (see Research Methodology). Eighteen 


months of recent demand figures were then collected and used 
2 





for model verification and forecasting comparisons. Various 
models were run with the sample of the FSC 6505 items selected, 
and the forecasts were compared to actual demands using the last 
two months as a primary test. The comparison of the forecasts 
to the actual demand were done in a managerial context; the 
assumption was made that the Supply Officer bought the amounts 
forecasted by the model, and the demand for the period is taken 
from the procured forecast. For each item, one of three con- 
ditions is possible; 1) The forecast will be greater than the 
actual demand, in which case there will be a surplus of material, 
and 2) The forecast will exactly equal the demand, in which case 
there will be no difference, and 3) The forecast will be less 
than the actual demand, in which case there will be a shortage 
of material. 

The evaluation of the forecasting model was done in two 
phases; an evaluation of the differences between the forecasts 
and actual demand, and an evaluation of the value of material 
associated with those differences. For each model, the total 
number of differences will be examined, with a view toward 
maximizing the total avallability (at the least cost of inventory). 
There will be two resulting "m M values, the value o£ 
material purchased that exceeds demand, and the value of material 
demanded that exceeds the forecast. The values are then 
summed to give the total material value. This results in 
weighing the value of material purchased that exceeds demand 
and the value of material demand that exceeds the forecast 


equally. 
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IV. RESEARCH METHODOLOGY 


DATA COLLECTION 

To test the stochastic models under consideration, a 
sample was required from the population of 1,006 items in 
FSC 6505 held by the Pharmacy at the Naval Regional Medical 
Center, Oakland. 

The demand information required for the line items was 
contained on the Retail Asset Stock Cards (RASC), compiled 
and sent to the Fleet Material Support Office, Mechanicsburg, 
Pennsylvania on a semi-annual basis. There is one card in 
the RASC deck for each item. Various information was con- 
tained in separate fields on each card as described in Appendix 
The FSC is contained in columns 8 through 20, Unit of 
mese in columns 23 and 24, 12 Month Demand in columns 57 
Arousa 62, and Unit Price in columns 70 through 76. The 
FSC 6505 section of the RASC deck submitted by the Regional 
Medical Center, Oakland in April of 1978 was reproduced by 
the Fleet Material Support Office and forwarded to the Post- 
graduate School for analysis. 

Two fields on the RASC cards were selected for analysis; 
the 12 Month Demand field and the Unit Price field. The 
total demand for the item over the twelve months preceeding 
the date of card submission was recorded in the 12 Month 
Demand field, and the Unit Price field held the system- 
designated unit price for the item. The demand figures were 


of primary interest, as the demand figures were used as inputs 
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to the models under analysis. 

To analyze the fields, one of the programs in the BMDP-77 
Biomedical Computer Programs (P-Series) was used. The BMDP 
series of programs, developed by the Health Sciences Computing 
Facility of the University of California, provides a set of 
easily understood and utilized computer programs that will do 
a wide variety of analyses, comparisons, regressions, and 
plots on data inputted by the user. The program used for the 
initial analysis was the P5D program, Histograms and Univariate 
Plots. The program's output is in the form of a histogram 
of the items in the field or fields designated, frequencies 
within the intervals of the histogram, percentages within 
the intervals, and the mean and standard deviation of the items. 
The 1,006 FSC 6505 cards were run with the required Control 
Language cards (Appendix 2), and the output of the program 
was obtained. 

The 12 Month Demand field was examined first. The mean 
of all 1,006 items was 548.009, and the standard deviation was 
1374.880. There were fifty intervals in the histogram, each 
interval containing 12 Month Demand Quantities. The interval 
boundaries were established by the P5D program, and all 12 
Month Demand quantities within those limits were contained 
in the interval. Twenty four intervals contained demand 
quantities. The first seven intervals accounted for 96.2$ 
of the demand (968 of the 1,006 items) with an upper limit 
of 2940 units of 12 Month Demand. Analysis of the first 
interval revealed that 34 of the 724 items in the interval 


held Ø demand in the 12 Month Demand field, and were 
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eliminated from the deck, leaving 690 items in the first 
interval. Since the remaining seventeen intervals (after the 
first seven intervals) contained only 38 items, or 3.8% of the 
total, a decision was made to disregard them in the selection 
of the sample. It was felt that the selection of an item 

from this group and its inclusion in the sample would not 
result in a representative sample; accordingly these items 
were treated as outliers. As a result, the sample would be 
drawn on a stratified sample basis from the first seven 
intervals. The interval limits, frequencies, and percentages 


are outlined below; 


12 Month Demand 


Intervals Frequencies Percentages 
per Lower 

420 0 690 690 VINEIS 73.8 

840 420 140 830 1: 9 881 
1260 840 38 86 8 4.0 22 
6 80 1260 29 Go 321 JO 
2100 1680 14 9 1.1 dieu 9733 
2520 2100 165 924 1.4 So 
2940 2520 10 934 1,3. 221006 0 


The Unit Price field was analyzed using the same program, 
P5D, Histograms and Univariate Plots. Again, the output was 
a histogram divided into fifty intervals, with each interval 
containing Unit Prices. Nineteen intervals contained Unit 


Price data, with 98.7% of the items contained in the first 
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ten intervals. Only twelve items, or 2.3% were located in 
the remaining nine intervals. The mean of the Unit Price 
data was $13.64, and the standard deviation was $27.05. 

The intervals, frequencies, and percentages of the first ten 
intervals are illustrated below: 


Unit Price 
* 


Intervals lec 15 Percentages 
Upper Lower 

10 0 689 689 69:3 0903 

20 10 125 814 125 915 

30 20 67 881 6.7 38.5 

40 30 36 217 36 o c 

50 40 28 945 2.28 94.9 

60 50 n 262 Ae. gene 

70 60 18 930 Im 98.4 

80 70 5 335 6 POM 

90 80 4 D = DES 

100 90 4 398 252 00.0 


* The range of the intervals are in $10 increments. 

As noted earlier, the 12 Month Demand field was selected 
as the source for the stratified sample, as the demand histo 
histories for the items would be used as the inputs into the 
inventory models. 

The RASC deck was sorted according to the upper and lower 


limits indicated on the 12 Month Demand histogram, resulting 
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in seven groups of cards, one group for each of the first seven 
intervals. Each interval card group was run through the BMDP 
P1D program of the BMDP Series, the P1D program being a Simple 
Data Description program. The program output was a complete 
listing of the group by item and number, actual 12 Month 
Demand, and Mean and Standard Deviation figures. Seven 
listings, one for each interval, resulted. 

The size of the sample was set at 100 items, approximately 
10% of the 968 remaining FSC 6505 items. The 100 items were 
selected from the seven intervals using the formula: 


X = 
7 X 100 = 2 


X = Number of items in the interval 


K 
N 


Total number of items in all seven intervals 


Z 


Number of sample items to be selected from the interval 
The problem then was to select items of the required sample 
Size from each interval. The procedure used was a modified 
random scheme using program ST1-04, Random Number Generator, 
Eme Of the Statistics Library programs sold for use with the 
Texas Instruments SR-52 Programmable Calculator. The program 
uses the Mean and Standard Deviation of the individual 
intervals to generate normal deviates within each interval. 
For instance, the required 75 random deviates for interval 0 
through 420 were generated from a normal distribution having 
mean 119.376 and standard deviation 109.567. As the normal 
deviates were generated, they were translated into individual 
sample items by matching with corresponding items listed on 


the print-outs from program BMDP P1D, Simple Data Description. 
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Each interval print-out had each item in the interval listed, 
along with the value of the 12 Month Demand associated with 
that item. As a normal deviate was generated, it was matched 
with the corresponding value in the 12 Month Demand field, and 
that item became one of the sample items. The required number 
of sample items were selected from each interval, and the 
result yielded a representative sample with the highest 
likelihood of selection near the midpoint of the interval. 

The items in each interval, number of sample items selected 


from each interval, and Mean and Standard Deviation of each 


interval are listed below; 


meer va | tee your Interval 

Upper Lower Std. Total Sample 

Limit Limit Mean Dev. Items Size | 
420 0 IS 6 690 75 
840 420 O TOC 140 15 
1260 840 dS CT GENE d 155 DIS 04 
1630 1250 1445.929 123.401 25 03 
2190 1680 JOSE 1526. 307 14 or 
2520 2100 2200300128050 ES onl 
2940 2520 2710000397784 10 On 
934 100 


After selection of the sample, 


the Means and Standard Deviations 


of the Intervals and the Samples were as follows: 





Interval 
upper 


420 
840 
1260 
1680 
2100 
2520 


2940 


atm C 


Lower 


Limit 


0 
420 
840 

1260 
1680 
2100 


2520 


Interval 


LAO 
592763 
10507708 
1445.929 
1853. 213 
2 69609 7 


2750.00 


Means 
Sample 
So 
3922 
993% 

1441. 


SE us 


2406 


2876 


94 


00 


25 


00 


00 


200 


00 


Interval 


Jo). 
Zee 
114. 
128 
20» 
T28 


Sa 


567 
183 
155 
401 
30 7 
850 


784 


* No Standard Deviation given as only one item in sample 


The Unit Price comparisons were as follows: 


ESC6505 Sample 
Items Items 
(1,006) (100) 
Mean Std. Dev. Mean Std Dev. 
2915.64 272085 $14.42 547.91 


The end result was a stratified sample numbering 100 items 
drawn from the FSC 6505 group numbering 1,006 items. 

Since the sample was identified by Stock Number for each 
the collection of demand history was considerably 


item, 


facilitated. At The Naval Regional Medical Center, the demand 
history is called the Usage History, and figures for January 

EU through June 1978 and July and August 1978 were collected. 
For each item in the sample, eighteen months of demand history 
two months of validation 


figures were assembled. In addition, 


figures were collected, (months nineteen and twenty) for use 


in forecast comparison. The eighteen months of history would 


be used to generate forecasts and establish demand patterns. 





Subsequently, the additional two months of history would be 
used to test the forecasts of demand for those two months. 
The eighteen months of demand data were transferred to 
machineable cards for use in the analysis. The Stock Number 
was included in Columns one through eight, while demand data 
were included in eighteen blocks of four columns each; all 
elghty columns of the card were used. The full deck of 100 
sample cards with the demand histories included bacame the 


Source Deck for computer analysis. 


ER MODELS 
Three models were selected for analysis; 1) A3RSR 
Method, 2) Exponential Smoothing Model with differing alpha 
values, and 3) Regression Forecasting. 
l. A3RSR 
The A3RSR Smoothing Method is actually an adaptation 
of a Tukey smoothing model into a forecasting model. A3RSR 
1s the name given to a smoothing technigue developed by 
rey and presented in the book Interactive Data Analysis by 
Donald R. McNeil.º As adapted for demand forecasting 
application, the smoothing is first applied to random data 
from a time series. The technique smooths the randomness 
into a relatively smooth curve, facilitating trend analysis 
and making forecasting possible. 
In actual operation, the smoothing is accomplished 
using running medians of three. The ith response in the 
series y(i-1), y(i), and Y(i-1) is replaced by the median 


of the three points. End points are treated by replacing 
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ABE by the median of y(1), y(2), and 3(y(2))-2(y(3)). The 
treatment of the last point in the series is handled in the 
same manner. It is this end point rule that makes the fore- 
casting possible. The use of median replacement in the time 
Bates data produces mesas, or pairs of adjacent points with 
a common value which is below or above the points on each side. 
The median replacement is repeated until the mesas are split 
and finally dissapear. The whole process can be viewed as 

a gradual erosion of the roughness of the original data in 
repeated small steps until no further smoothing by medians 
is possible. In fact, the '3RSR' part of the name refers to 
"medians of 3, repeated until convergence, split, repeated 
mail convergence".? The "A" part of the name refers to the 
use of medians as a basic component of the smoothing as 
opposed to a "B" technique, explained by Mr. McNeil in his 
book. 8 

At the Naval Postgraduate School, the A3RSR program 
is part of the OA366f@ APL Public Library on file in the IBM 
360/70 Computer in the W.R. Church Computer Center. Data 
were transformed into the correct APL format for the program, 
the program run, and the forecast: gathered as follows; 

The Source Deck (demand history cards) was read into 
the computer, and set up by a FORTRAN program into a file named 
FTO4F001. The APL mode was entered, and the file FTO4F001 
was read and transformed into an APL formatted file named 
FTO8FO01. File FTOS8F001 was renamed 'DEMAND', and the Stock 
Number deleted from the file. The final form of the file 


named 'DEMAND' was a line of demand data for each item in 





an APL format, where each item is accessed by requesting the 
file name, then the numerical number in the sequence (the 
third item is requested by the command DEMAND (3;) ). The 
Public Library OA3660 was copied into the APL workspace, and 
each item in the 'DEMAND' file was run through the A3RSR program. 
The result was assigned the value of 'C' (e.g. Ce-ABRSR DEMAND 
(21;) ) and 'C' was printed out. The last numerical value in 
the print-out of 'C' was used as the A3RSR forecast for the 
next period (Appendix 3). 


2. Exponential Smoothing 


The Exponential Smoothing technique, set forth in a 


ERR entitled Statistical Forecasting for Inventory Control 
by Robert G. Brown, published in 1959, is essentially a 


9 It uses the 


weighted moving average method of smoothing. 
Semana tor the present period (d+), and the forecast for the 
last period uU and a constant labled '«' to develop 
poaorecast for the upcoming period (UL). Mathematically, 


the formula appears as follows; 


Ue = (4) (do) + (1-4) (Up) + 
U, = Forecast for upcoming period 
as = Demand for the present period 
Ur-1 = Forecast for the last period 


As the new demand figures for the present period (d+) 
become known, they are worked through the formula with the 
PEN canta and the forecast for the last period (U.3) to 
Pem@etate the new Forecast (U_). Exponential Smoothing is 


especially adaptable to data processing in that only the value 


2) @ 





S fl; S Or St (Urey) need be stored for each item, 
making for a minimum amount of storage.tl 

One facet of the Exponential Smoothing model is the 
fact that the value for « can be changed to give different 
weights to the demand in the present period (d,) and the 
me sr tor the last period (Uf-1]). Analysis of the present 
model in use at the Naval Regional Medical Center, Oakland 
reveals that the model is in actuality an Exponential Smoothing 
model with an « value equal to 0.333. See below; 


BE 2T1OPT1) 
= OP2 is equivalent to (CD) + -$-(0P1) OP» 


3 
Sua?) = Cure Demand SSME period) 
OC = Forecast (order Point) Last Period 
OP, (UL) = Forecast (Order Point) for upcoming period 


This ability to change the value of « raises the 
question of the possible benefits of changing the values and 
thereby changing the weights. In the course of the study, the 
decision was made to analyze the effect not only of the current 
« value (0.333), but also the other & values. The values 
A ere 07001, 0.01, 0.025, 0.05, 0.1, 0.2, 0.333 
(present system), 0.5, and 0.7. The values lower than 0.5 
put ever increasing weight on the forecast for the last period 
(Ut-1), the value of the 0.5 puts equal weight on both the 
AMEN demand (d+) and the forecast for the last period (Ur-]), 
and values greater than 0.5 put the preponderance of weight 
EN ucocurrent demand (dec). By using these values of 


in the analysis, the effect of changing weights in the formula 
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can be assessed, and the effect of changes in the present 
system formula can be illustrated for comparison with the 
other models. 

In actual operation, a FORTRAN program was written to 
develop the forecast for the nineteenth and twentieth months 
(Appendix 4). Once the forecasts were printed by the program, 
they were compared to the actual demands for the periods, and 
the deviations were recorded. 

3. Regression Forecasting 

Regression forecasting employs standard and well-known 
regression formulas with time-series analysis.-? Standardized 
formulas use the demand data inputs to give the Slope and 
Y-Intercept of the data set, while the time-series element is 
brought in by including time data (y(i).b The regression 
forecasting rule used in this study is based on the format; 

y(i+l) = Y-Intercept + Slope (y(i))!# 
The formula is derived from the standard linear regression 
model in which a straight line is fitted through data points 
using the method of least squares. The data thereby yield 
estimates of the Y-Intercept and Slope which are then used 
in the actual forecast. 

To effect this procedure, the eighteen month demand 
data were processed through the 'LINE' program of the OA366f 
APL Public Library, which resulted in Slope and Y-Intercept 
figures for the demand data. The APL commands to use the 


'LINE" program are, for example; 
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Dp DEMAND (16;) 
W K LINE D 
SLOPE: 1.052 s TIP RCEPT:. 104932 

The first command assigns to 'D' the values of the eighteen 
demands in the 16th line of the 'DEMAND' file, while the 
second command runs the values of 'D' through the 'LINE' 
program and asSigns the estimated values of the Slope and Y- 
Intercept to 'W', which is then printed out in the form of 
line 3. The forecasts for both the nineteenth and twentieth 


months are developed in this way, and are compared to the 


actual demands for these periods. 


Gee ANALYSIS 

As stated earlier, the eighteen months of demand history 
for each sample item were entered on machineable cards with 
Ee Stock Number for each item, which resulted in a 100 card 
Source Deck to be used as the input into the various stochastic 
models. 

To get a pictorial view of the eighteen month demand for 
each item, another program in the 0A3669 APL Public Library 
was used, the 'SCAT' program. This program plots a graph of 
the data points, placing the demand values on the Y axis and 
the individual months on the X axis. A demand scattergram 
is created by the program that visually depicts the demand 
values over the months included. Each item in the sample 
was run through the 'SCAT' program so that an idea of the 
overall pattern of the demand points for each item could be 
seen. 


The A3RSR forecasts were developed as explained in the 
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A3RSR model section for each of the two months for which data 
was available. The forecasts were listed along with actual 
demand recorded for each month, and notations were made of 
the differences. 

The Exponential Smoothing forecasts were developed for 
each value of A from the FORTRAN program as outlined in 
Appendix 4. The forecasts were listed along with the actual 
demands recorded for each month, and notation made of the 
differences. 

The Regression Forecasts were developed as outlined in 
the Regression Forecast section for each of the two months 
for which data was available. The forecasts were listed along 
with the actual demands recorded for each month, and notation 
made of the differences. 

Differences between forecasts and actual demand for each 


item were illustrated by the following examples; 


Brock No. Forecast Demand Difference 
B1900 35 30 +5 
Zm 182 1 29 27 00 
2258100 73 Je -25 


The differences are explained in a managerial context: 
The assumption is made here that the Supply Officer purchases 
the amount indicated by the forecast to meet the upcoming 
demand. If the purchased (forecasted) amount exceeds the 
demand (as in the first case above), a surplus exists at the 
end of the period. If the purchased (forecasted) amount 
exactly meets the demand (as in the second case above), 


there is zero difference. If the purchased (forecasted) 
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amount is less than the demand (as in the third case above), 
a shortage exists at the end of the period. For each model 
evaluated, the number of surpluses, zero differences, and 
shortages will sum to 100, the size of the sample evaluated. 
The measure of differences forms the basis of the 
quantitative cost evaluations for each modei. The differences 
for each item are multiplied by the Unit Price of that item, 
giving either the value of material that would have been 
required to meet the demand if the forecast had been accurate 
or a value of material demanded that exceeds the forecast. 
If a model were 100% correct in forecasting demand, there 
would be no dollar value differences associated with the 
model (dolar value differences as defined in the preceeding 
sentence, of course). The only dollar values would be the 
value of material purchased that exactly met the demand. It 
should be recognized that no model will be 100% correct (due 
to the stochastic nature of the demand), so there will 
inevitably be values of material purchased that exceed demand 
and values of material purchased that exceeded the forecasts 


associated with each model. 


V. FINDINGS AND CONCLUSIONS 


As stated earlier, forecasts were developed from each 
model using the formulas described in the Research Methodology 
Section. If a forecast was in decimal form, it was rounded 
to the next whole number (e.g. 1.5 was rounded to 2.0); this 


was done to conform to the managerial context of the analysis. 
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The manager, under this assumption, purchases the forecast, 
and rounding the decimal forecast is required to determine 
the number of units to procure. 

The forecasts were compared to the actual demands for 
the two months under study, and the differences enumerated 


for each model. The differences are listed below; 





Ich 20th Per- 

Model Month Month Total Cent 
A3RSR 

(+) 61 42 103 a 

(9) 12 LS 25 12055 

E) ZEN 45 => Bom) 

100 100 200 100.0 


sm. (0.001) 





(+) 58 48 106 53.0 
(9) 04 10 14 07.0 
i) 38 42 O 40.0 

100 100 200 100.0 


ER sm. (0.01) 


(+) 60 46 106 S 
(2) 02 06 08 04.0 
E 38 48 86 43.0 

100 100 200 190070 


B Sm. (0.025) 





(4) 62 51 113 56.5 
(9) 03 02 05 32 55 
(>) 35 47 _82 41.0 

100 100 200 200,0 
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Model 
Ex.Sm. 
(+) 
(9) 
e) 


E Sm. 


(+) 
(9) 


Ec Sm. 


(+) 
(9) 
©) 


coos, Sm. 


cn) 
(9) 
(=) 


Ex. Sm. 


(+) 
(9) 
>) 


19th 
Month 


(0.05) 


63 

05 

p 

100 
40151 


67 


(0.333) 
73 
02 
Bs 
100 
(95) 
75 


00 


25 


100 


20thn 
Month 


49 
02 
49 


100 


47 


48 
07 
45 


100 


52 


30 


Total 


Jio 


07 


124 


200 


Per- 
Cent 


506: 0 
OE 


40.5 





100.0 


2400 
USO 
30.5 


10020 











00.0 





19tn 20th Per- 
Model Month Month To cal Cent 


E Sm. (0.7) 





(+) 72 41 113 56.5 
(9) 01 13 14 07.0 
(=) 27 E ES 36.5 

100 100 200 100.0 


Regress. Fore. 


(+) 63 35 101 5075 
(9) 08 10 18 09.0 
(=) _29 _52 _81 _40.5 

100 100 200 100.0 


On the basis of the absolute values of the total positive 
and zero differences (which reflect the ability of the system 
to fill the demand when first requested), the Exponential 
Smoothing model with an alpha value of 0.2 forecasted best 
both within the Exponential Smoothing group and against the 
other two models. Viewed in the context of the restricted 
number of months compared (two), the differences in the 
models were not considered to be alarming, in that the total 
Variance range between the two models was only 6.5%. If 
more periods had been avallable, it would have been reasonable 
to expect that the percentages would change, and a better 
evaluation of the first pass availability would be possible. 

To assess the differences over a longer period of time, 

a Smaller sample consisting of the mean items from each 
interval was run with the A3RSR, Exponential Smoothing 


(alpha= 0.2), and Exponential Smoothing (alpha= 0.333) models 


BI 








per 
rev Vam HI 


Bea period of sixteen months. The full twenty months of 
data was used, but due to a pecularity in the A3RSR program, 
the first four months of demand data could only be used to 
develop forecasts. The running of the entire 100 items in 
the sample for a period of sixteen months was considered to 
be a task of too great a magnitude given the time constraints 
of the study, so the mean items in each interval were chosen 
as a representative sample of the 100 items. The results 


are illustrated below; 


A3RSR E MI) BA Sm 0922) 
ED 55 5 (5) 
(9) 04 (g) 02 (g) 01 
SS (a ER SE 
12 112 IBI 


AS can easily be seen, the differences are consistent 
with the results on the larger sample. Model performance 
over time appeared to be not greatly different from per- 
formance over the two month period. 

The dollar values derived from multiplying the differences 


by the Unit Price of the items were summed and are presented 


below; 
IER Ai 
Models Month Month Total 
A3RSR 
E) 2020/2. 60 $4,410.74 sL1,437.94 
e) 2252.32 223290522 6,486.14 
Rota). 297 219.252 $8,643.96 S77 923.48 
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Models 


Er Sm. 


i 
ir) 


Metal 


EX. Sm. 


(+) 


Total 


Metal 


E Sm. 


G) 


Total 


Co On!) 


uo) 


05025) 


MOTOS 


W) 


m2) 


19th 
Month 


5575555 dl 


3,463.86 


IOS 07 


220 30 


20 D32 


Saal 90 822 


o o l7 


2 SB. AA 


Siero oO. Lo 


OO US 91. 


25572. 15 


Sooo 09 


09 23102065 


ino 1-753 


SS 


096.45 


MVS O17 


555 265. 02 
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Month 


5 999102 


37030.41 


21200629223 


247 2660.25 


MIOS O 


Sur 39. IL 


san 321256 


5 12.71 


5106058343. 27 


55,120. 74 


O Le 


2072922995 


95,184.84 


EE 049.9 


$8,949.39 


S57 P3730 


o nO E 


IS 
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Total 


S10 7637223 


_ 9,094.27 


00 


S EO 


87994.15 


SS 


1027795275 


MIES 


s 159 700153 


Su pas» 


6,724.40 


S75 3505 


Su 116.49 


2460426 


ni 65 0062. 77 


2.093275 


VÁ SS: 


6853435 





Models 


BRs Sm. 


(+) 
(=) 


Total 


cece Sm. 


i 


Total 


Iu Sm. 


i) 
(isi) 


Metal 


Regress. 


(+) 
(-) 


Total 


o) 


0:5) 


(02 7) 


Fore. 


Month 


57,430. 


mer 


3,329. 


SOL 


28,924. 


ZO « 


06: 527. 


58,286 


78 


30 


1320: 


08 


48 


w229: 


45 


45 


47 


90 


BD 


5 


70 


42 


IEEE 


¿L 2 


Month 


5,021. 


23,493 


$4,696. 


38,035. 


$4,549 


3287339 


230857 


Sonos 


22 


40 


an 12, 


262 


46 


DO Oe 


90 


36 


226 


no 


4,290. 


„ehe 


220 


86 


47173: 


83 


Total 


12,452. 


$17,244. 


ous 7225 


$18,354. 


Susp en 


S poe 


SA 


056/025 


00 


70 


au 192% 


70 


„94 


DL o 


87 


61 


73 


62 


5,642. 


33 


67 


"357. 


28 


95 


In terms of absolute dollar value, again the Exponential 


Smoothing model with an alpha value of 0.2 performs the 


best (gives the lowest dollar values) 


other other models. 


the value differences, 


in relation to the 


When evaluated as to the magnitude of 


the models perform nearly equally. 


Changing the values of alpha in the Exponential Smoothing 


models produced differing dollar values, which proved to 


be lowest where the alpha value equalled 0.2. 
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These 


202509 


79350 


18450 


27550 


16650 


empirical results can be best seen 


in the graph below 
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An overall graph of the Exponential Smoothing model with the 





different values of alphas and other stochastic models proved 


to be very interesting; 


20250 = 

i 

E eco One) 

| (02015 

| O 

I 
19350 -- (0.7) 

| o 

| 

| 

N (0.025) 

| O 
18450 = SS) 

| 

| 

o A3RSR CON OS) 

I O O 

| REGRESS 
M0 |. OZE Ob 

I 

| 

| o (9.333) 

N OSO.) 

O O 
teca 

0.0 25 5 Ü ge 10.0 12.5 
Other Models Total Dollar Values 

A3RSR S 923018 
Reg. Fore. Su 002995 


AS can easily be seen, not all the Exponential Smoothing 
variations did better than the A3RSR and Regression Forecasing. 
Only certain values of alpha produced. total dollar values in 
the Exponential Smoothing model better (Lower) than both the 
other models (0.1, 0.2, 0.333). Even for those values below 
the values of the A3RSR and Regression Forecasting models, 

the differences were not large enough to give unqualified 


approval to the Exponential Smoothing model above all others. 
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One notable fact is that the total value for the model currently 
in use at the Naval Regional Medical Center, Oakland (Expo- 
ential Smoothing with an alpha value of 0.333) is lower in 
total dollar value than the A3RSR and Regression Forecasting 
models, and falls below both of them on the graph above. 

The Value Difference/Ex. Sm. (0.2) Value of the models 
are presented below; 

A3RSR Roque ssicalnorecasE 


6.33 4,4 


oo 


Here, as in the differences comparison, the percentage range of 
variation was within 6.5%, leading to the conclusion that 
there is no real great difference between the three models in 
their performance. 

The valuation of the differences (total dollar value) 
found in the mean items drawn from the sample and run with 
the forecasting models tended to support the results discovered 
in the 100 item sample. The period was sixteen months, but 
the results followed the same pattern as the two month 


primary test; 


A3RSR BDC: OMV 8) Es memo (02) 

= 59,921.31 (+) $9,556.03 2221-9002 502 

() 9,645.91 (9 ‚238289 ee 
Mercal 519,567.22 MS OZ 180 214,.01 


Here again, in absolute dollar value terms, the Exponential 
Smoothing model with an alpha equal to 0.2 proved to be the 
best. The magnitude of the differences, however, was not 


large enough to give an overall endorsement to the Exponential 
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Smoothing model. 

The finding that there was no large difference between 
the performance of the three models led to an attempt to 
discover why. A more significant variation in model performance 
had been expected due to the difference in forecasting approaches 
inherent in each model. Yet it appears that no matter what 
model was used, the results were nearly the same. 

In an effort to explain the lack of difference in model 
performance, the patterns of demand for the interval items 
were examined. Previously, the eighteen months of demand had 
been graphed using the 'SCAT' command and a line fitted to 
the points using the 'LINE' command. Residual points were 
plotted around the fitted line by using the 'SCAT' command 
after fitting a line to the random demand points. The 
result was a plot of the points around the fitted line, 
which graphically illustrated the degree of dispersion in 
the data. 

The spread of points around the fitted line proved to 
be quite wide, which is to say that the data in general have 
a large inherent variance. Large variance in the demand data, 
greater than what would be normally expected, may very well 
explain the robustness of the different models. For best 
performance, the Regression model should have data that 
conform well to a linear structure. The closer the data 
point fit to a trend line, the better the model performance 
will be. The Exponential Smoothing model is set up to nandle 
more variable data; i.e., is more adaptive in nature. The 


A3RSR technique attempts to break down outliers and wild 
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variation to discover an underlying structure. Beyond a 
certain point, however, the amount of variability and/or 

bias may be too large as to completely obscure any reasonable 
attempt to discover underiying structure. 

Further investigation was thus undertaken to try to 
Understand better all the factors impacting on the system. The 
demand was originated by the Pharmacy when they ordered the 
FSC 6505 items from the Supply Department, so the origins of 
the orders inside the Pharmacy were investigated. 

The Pharmacy Officer was interviewed concerning the 
ordering procedures in effect within his Division. Supple- 
mentary information was also sought concerning the dollar 
value and unit amount of material surveyed (destroyed) on 
a monthly basis. The reason for the survey information was 
to discover if the excessive fluctuation observed in the data 
had caused the Pharmacy to overstock material, which then 
expired and had to be destroyed. If the Pharmacy overstocked 
material due to excessive random demand, then the demand 
histories used by the Supply Department would cause a cor- 
responding overstockage in the Supply Department. The Supply 
Department surveyed $1,042.28 worth of expired material from 
January through June of 1978, so it was hoped that a comparison 
could be made to see if overstocking due to excessive demand 
affected both areas. 

Within the Pharmacy, the ordering of FSC 6505 items is 
dependent on two factors; 1) The demand for material at the 


Pharmacy window (from 'customers' of the Pharmacy), and 2) 
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The determination by the person doing the ordering for the 
Pharmacy of the requirements for material to support the 
demand at the window. The demand excountered at the window 
is certainly the central source of randomness in the system; 
everything else should be dependent upon it. Even though 

the window demand is random, a means of measuring that demand 
should exist to give the person doing the ordering some idea 
of the demand history of that item. 

Unfortunately, no means currently exists for measuring 
the window demand within the Pharmacy. No system is in 
effect to give any idea of the demand history for an item. 
The reason for this (given by the Pharmacy Officer) is that 
only a week supply of FSC 6505 items are kept on the shelf 
in the Pharmacy, so a system for recording demand is not felt 
necessary. A system has been proposed whereby the Data 
Processing Division will provide to the Pharmacy on a line- 
item basis a four month demand history. Within the Pharmacy, 
it was felt that this report would not be used even 1f 
provided. The Pharmacy personnel felt that they were 
performing at an optimal level, and did not need a demand 
history report. 

The process of ordering from the Supply Department was 
then investigated. The ordering system is subjective in 
nature; the person doing the ordering keeps a card on each 
line item used, and records on that card the date and quantity 
to be ordered that he subjectively determines he needs. 

Each time he orders, he consults his previous order amount, 


me orders what he thinks will fill the demand within the 
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next week. The person doing the ordering has no idea what the 
real past demand has been at the window, so his order amounts 
are truely subjective and tend therefore to be highly biased. 
The records concerning material surveyed within the Pharmacy 
are kept for only narcotic items, as this is required by 
Federal law. Other FSC 6505 material surveyed is either 
thrown away without recording, or, in the case of material 
purchased from commercial sources, is given back to the drug 
Salesmen for replacement. No estimation of the dollar amounts 
or volume of surveys Can be made, as no records are kept. The 
Pharmacy has no objective system to determine how great the 
demand is at the window, how great the demand has been in the 
past, or how much material is surveyed each month. 

The lack of data collection within the Pharmacy concerning 
demand, the lack of survey information, and the subjective 
Ordering system certainly contribute to excessive fluctuation 
in the demand data, with the effect of obscuring the true 
demand pattern. Three factors are needed in the Pharmacy to 
establish an objective ordering system; 

1) Measurement of the random demand encountered at the 

Pharmacy window 

2) Measurement of the amount of material surveyed 

whether required by law or not 

3) Measurement of the amount of material needed from 

the Supply Department to meet the demand encountered 
at the window 


Within the Pharmacy, there is no objective system to satisfy 
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the requirements for these measurements. Therefore the 

amount of material ordered from the Supply Department can 

vary widely on a number of factors. Seasonal demand may 

cause temporary stockouts in the Pharmacy that will cause 

the person ordering material to overstock, different personnel 
may be doing the ordering, excess material ordered one month 
may be depleted over a number of months distorting the demand 


figures, and the desire to conform to funding limits may 





influence order amounts. In a system such as this, it is 
easy to understand why no model can perform appreciably 


better than any other. 


VI. RECOMMENDATIONS 


Since excessively variable data renders the performance 
of differing forecasting models nearly equal, steps should 
| be taken to reduce the variation in the data, or at least 
| account foxr it better. The demand data from the Pharmacy is 
| dependent on three factors; 
1) The random demand encountered at the window 
2) The amount of material surveyed 
3) An objective system to determine requirements for the 
upcoming period 
The Pharmacy should establish a system to measure the 
demand for drug items at the window. This is a first step 
in establishing an objective system for the ordering of 
material. The building of a demand history will both aid in 


ordering material and aid the Pharmacy personnel in 
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understanding the material they handle. With a good demand 
history they can identify seasonal items, items with changing 
demand, and items needing special handling. 

The amount of material surveyed should be recorded 
within the Pharmacy for their internal use. Survey data, when 
used in conjunction with the demand encountered at the window, 
gives a better picture of material requirements to the personnel 
doing the ordering for the Pharmacy. Secondarily, accurate 
s y records aid in drug control, reducing the possibility 
that theft may be explained as material surveyed. 

To begin establishing an objective ordering system, the 
Pharmacy should work with the Data Processing Division in 
setting up the demand history program that currently is being 
considered. The system should give the Pharmacy personnel 
average demand data for the past four months, which will at 
least give the personnel doing the ordering an idea of how 
mcuh material should be ordered. The user involvement of 
the Pharmacy in the program development will help in assuring 
that a viable program will be developed, one that will be 
utilized rather than left on a desk. Development of this 
program will be an important step in reducing demand data 
variance, and is best undertaken by both departments involved. 

Average demand data is more objective than an estimation 
of requirements based on previous order amounts. If average 
demand appears to give enough information for ordering 
purposes, no further system development need be undertaxen. 


If excessive stockouts or surpluses result, further refinements 
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can be undertaken to refine the ordering system to make it 
more responsive. Whatever system is implemented, it should 
be at least as good or better than the subjective system in 
operation at the present time. 

Once an objective system is in place and the randomness 
in demand data reduced, the study should be rerun to evaluate 
the models under the new demand circumstances. With less 
random and more representative demand data, model performance 
differences should be more pronounced, and a better evaluation 
of the best model to use for demand forecasting should be 


possible. 
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APBENDIX A 


RASC CARD FORMAT 


Document Identifier (BA5) 

Condition Code 

VOSL Item Leadtime in Months 

Overflow Indicator 

National Stock Number (NSN) 

Logistic Management Code "F" or "G" (9M Cog) 
Blank 

Unit of Issue 

On-Hand Quantity 

Due Quantity 

Acquisition Advice Code 

Back Order Quantity 

* To denote gas cylinder item 

Cylinder Pool Requirement/VOSL Reorder Point 
PWRS Reservation Quantity 

VOSL QMC (One Decimal) 

Necoumte/ Cosi ıB MIW, 9A, 9C, etc.) 

Demand (12 Months) 


Demand Frequency (Each Quarter) 


Roüting rdentifler (Erom) 
Unit Price 


NSO Quantity 
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APPENDIX B 


PROGRAM P5D, HISTOGRAMS AND UNIVARIATE PLOTS 


fee control Language 


mee identification Card 
T // EXEC BIMED, PROG=BMDP5D 


ERR SYSIN DD * 


program PSD Control Language 


4. /PROBLEM Title is 'FSC 6505 Line Item Histogram Plot'. 


Er /INPUT Variables are 2. 

6. Formas (56xXx6 0 7 F/.2,4X)'". 
7. /VAR Names are '12MODMD' ‚Uprice. 

oe Blanks are missing. 

9.  /PLOT Tyve=Hist, Norm, Cum, Chist. 

TO. Variable = '12MODMD', Uprice. 

IM /END. 


Data Cards 
Deck of 1,006 Line Item Cards 
End Card 


ee / / 
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APPENDIX C 


A3RSR FORECAST PROGRAM 
This program uses the 100 card Source Deck as input. The 
FORTRAN Program converts the Source Deck into a file named 
FTO8F001, which is then converted into APL format by the 
APL commands. A file named 'DEMAND' is created by the APL 
commands, and the individual items in the file are run 
through the A3RSR program of the OA3669 Public Librarv. 
l. Read Source Deck into outside card reader. (name 
mente FTO4F001') 
2. Run source Deck through following FORTRAN program: 
REAL DATA (100,22) 
N=100 
READ(4,1) ((DATA(I,J), J=1,22), I=1,N) 
IESEORMAT(IX,E2.0,1X,19E4,0/72E4,.0) 
DO 10 I=1,N 
10 WRITE (8,2) (DATA(I,J) ,J=1,21) 
ERG EMAT (7 3.0,1X,25.0, 198520) 
SLOP 
END 
3. The APL mode is entered, and the following APL program 
created; 
IC OPY EC MSTO CMSRD 
DEMAND l CMSREAD 'FILE FTO8F001' 


STOCKNO Emo 2 DEMAND 
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) WSID INVENTORY 

) COPY 20A3660 ESSENTIALS 
4. The APL Program reads the FORTRAN file FTO0S8F001 and 
converts it to an APL file named DEMAND. 
5. Individual line items are selected from the ' DEMAND' 
file and run through the A3RSR program; 

Di €- DEMAND (1;) 

C «— ABRSR Dl 


€ 


6. At the command 'C', the A3RSR forecast is printed out, 
and the last value is taken as the forecast for the coming 


period. 
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APPENDIX D 


EXPONENTIAL SMOOTHING PROGRAM 


This 
the 
use 
ie 
and 
fore 


Job 


program developes Exponential Smoothing forecasts using 
Source Deck as an input. Alpha values are adjusted by 
of the command 'Al=X' where X is the alpha value desired. 
program develops forecasts for two months at a time, 

uses the means of the first three demand values as the 


east for the last period in the first forecast calculation. 


Control Language 


Ike 
z 
3: 


4. 


Identification Card 
LU EXEC FORTCLG 
7/ GO. SYSIN DD * 


A 


Exponential Smoothing FORTRAN Program 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 


2911 


DIMENSION x020 739 990 STOCK (100), 2(20,100) 
REAL*8 STOCK 
DO 10 I=1,18 
X(I)zI 
10 CONTINUE 
pon 7=1,100 
Ber 05,100 ,END=-INESEOER II) ,(Y(K,J),K=1,20) 
100 FORMAT(A8,18F4.0/2F4.0) 
20 CONTINUE 
99 CONTINUE 


WRITE (6,610) 





WEZ 610 FORMAT(' STOCK NUMBER PREDICTION') 


0013 powS9OLI-l,100 

0014 Al=.1 

0015 Al=1.-Al 

0016 Gee AE M A2 AY (LE) + Y (2,1) +Y(3,1))/3.) 
0017 DO 300 J=2,18 

0018 Z(J,I)=A1*Y(J,I)+A2*Z(J-1l,¿ć{I) 


ES 300 CONTINUE 
0020 PE O0 LOC (ZO), AL, II 


ERRA GOO FORMAT (2X,A8,8X,F8.1,10X,F5.0,I5) 


0022 DNE 

0023 A2=1.-Al 

0024 Z O U DO a) (2,1) +1 (8 ,1))/3.) 
0025 bemdamnegs» 8 

0026 A AO SEE) 


0027 400 CONTINUE 


0028 Ene o do) AS TOCA z(18/1), AI;II 

0029 DNE M 

0030 Por 

0031 A EY (2 1) +1 (3,1))/3.) 
0032 DO 310 J=2,19 

0033 POUR TV (3,1) A222 (3=1 1) 


MOBA CONTINUE 


0035 rus (o 600) STOCK (L) Z(197I); AL, ,II 

0036 Al=.05 

0037 A2=1.-Al 

0038 AAA, IEA A (SL, LAY (2,1) +1(8,1))/3.) 
3039 DO 410 J=2,19 
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0040 


0041 


0042 


0043 


0044 


0045 


410 


200 


DADA LAA AS AS LL) 

CONTINUE 

Whiten 500) STOCK), 2(19,1), Al, LI 
CONTINUE 

RETURN 


END 


ll 





IE. 


IL. 


i. 


I. 


14. 
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